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The  effects  of  absorbed  moisture  and  temperature  on 
the  complex  moduli  of  composite  laminates  are  investigated 
and  the  mechanisms  of  moisture  diffusion  in  a  lamina  are 
also  analyzed . 

First,  the  variation  of  the  complex  moduli  of  epoxy 
in  terms  of  temperature  and  moisture  content  are 
experimentally  determined.  Then,  the  hygrothermal  effects 
on  the  complex  moduli  of  composites  are  derived  by  using 
the  complex  moduli  of  the  matrix,  mi cr omechanical  formulas 
and  experimental  data.  Only  the  hygrothermal  effects  on  the 
complex  moduli  of  pure  epoxy  need  to  be  experimentally 
determined  since  these  effects  on  the  fibers'  properties 
are  negl igib 1 e . 
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In  addition,  the  effects  of  hygro therma 1  environment 
on  the  stress  field  and  material  damping  of  general 
laminated  composite  plates  are  analyzed.  It  is  shown  that 
hygrothermal  stresses  induced  directly  by  moisture  and 
indirectly  by  material  property  changes  can  be  very  high, 
but  the  effects  on  damping  are  less  pronounced. 
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CHAPTER  1 
INTRODUCTION 


1  .  1  General  Introduction 


The  introduction  of  advanced  composites  in  aerospace 
applications  has  led  to  an  extensive  study  of  their 
mechanical  behavior.  The  amount  of  experimental  and 
theoretical  findings  of  composite  material  researchers  made 
during  the  1960's  was  so  vast  that  Broutman  and  Krock  [1] 
needed  eight  volumes  to  edit  a  summary  of  the  resulting 
knowl edge . 

The  interest  in  composite  materials  arose  from  their 
ideal  characteristics  for  aerospace  structures.  Replacement 
of  the  commonly  used  aircraft  material,  aluminum,  by  high 
strength/density  ratio  and  versatile  composites  can  lead  to 
a  theoretical  60%  weight  reduction  [2,  p.  22].  Due  to  such 
benefits,  lower  costs  and  better  understanding  of  their 
mechanisms,  the  use  of  composite  materials  has  been 
increasing  slowly  but  steadily. 

Exposure  of  aircraft  structures  to  high  temperature 
and  humidity  in  the  environment  and  the  tendency  of 
composites  to  absorb  moisture  gave  rise  to  concern  about 
their   performances   under   adverse   operating   conditions. 


Therefore,  considerable  work  has  been  done  to  understand 
the  effects  of  hygrothermal  environment  on  the  mechanical 
behavior  of  composite  materials. 

1  .  2  Moisture  Diffusion 

In  a  1967  study  on  the  effects  of  water  on  glass 
reinforced  composites.  Fried  hypothesized  that  water  can 
penetrate  the  resin  phase  by  two  general  processes,  by 
diffusion  through  the  resin  and  by  capillary  or  Poiseuille 
type  of  flow  through  cracks  and  pinholes  [3].  But  no 
mathematical  theory  was  presented.  Later,  investigators 
established  that  the  primary  mechanism  for  the  transfer  of 
moisture  through  composites  is  a  diffusion  process  and 
adapted  the  general  theory  of  mass  diffusion  in  a  solid 
medium  to  moisture  diffusion  in  composite  materials.  The 
transfer  of  moisture  through  cracks  is  a  secondary 
effect  [4,  5].  Experimental  data  indicate  that  for  most 
composite  materials,  the  diffusion  of  moisture  can  be 
adequately  described  by  a  concentration  dependent  form  of 
Fick's  law  [4-10]. 

1 . 3  Hygrothermal  Effects 


The   degradation   of   mechanical   properties   of   glass 
reinforced   plastics   exposed   to   water   has   long   been 


recognized  by  marine  engineers  who  use  "wet"  strengths  in 
the  design  of  naval  structures  [3].  Requirements  in 
aircraft  structures  are  more  stringent.  The  mechanical 
properties  of  materials  used  in  aerospace  applications  must 
be  completely  characterized.  Therefore,  the  effects  of 
hygrothermal  environment  on  the  elastic,  dynamic,  and 
viscoelastic  responses  of  composites  have  been  studied.  To 
date,  the  effects  of  moisture  and/or  temperature  on  the 
following  performances  have  been  investigated:  tensile 
strength,  shear  strength,  elastic  moduli  [3,  11-14], 
fatigue  behavior  [15-17],  creep,  relaxation,  viscoelastic 
responses  [18-20],  dimensional  changes  [21],  dynamic 
behavior  [22],  glass  transition  temperature  [23],  etc. 

Only  tensile  and  shear  strengths  and  elastic  moduli 
have  been  thoroughly  studied  by  many  researchers.  But  data 
on  the  other  properties  are  more  limited  and  hence 
inadequate  to  constitute  a  good  design  data-base. 

1 . 4  Scope  and  Methodology 

The  present  investigation  is  a  combined  theoretical 
and  experimental  work  and  is  concerned  with  predicting  the 
hygrothermal  effects  (below  the  glass  transition 
temperature)  on  the  complex  moduli  of  composite  materials. 

This  program  is  undertaken  by  carrying  out  the 
foil  owing  s  teps : 


i)  The  complex  moduli  of  epoxy  matrix  in  terms  of 
temperature   and   moisture   concentration   are   obtained   by# 
using  experimental  tests  and  theoretical  expressions. 

ii)  The  effects  of  temperature  and  moisture  on  the 
complex  moduli  of  unidirectional  composites  can  be  derived 
by  using  the  complex  moduli  of  the  matrix,  micromechani cs 
formulas,  and  experimental  observations.  In  addition,  we 
neglect  the  hygrothermal  effects  on  the  fibers. 

iii)  The  effects  of  hygrothermal  conditions  on  the 
stress  field  and  the  material  damping  of  some  general 
laminated  composite  plates  undergoing  simple  hygrothermal 
loadings  are  analysed. 

1 . 5  Dissertation  Lay-Out 


Right  after  the  introduction,  the  mechanism  of 
moisture  diffusion  is  described  in  Chapter  2,  where  the 
absorption  of  moisture  through  thin  composite  laminas  is 
analyzed  in  detail. 

The  complex  moduli  of  unidirectional  composites  are 
defined  in  Chapter  3.  Sections  3.3  and  3.4  give  the  microm- 
echanics  formulations  of  the  elastic  and  complex  moduli  in 
terms  of  the  constituent  material  properties.  The  damping 
of  composites  based  on  the  dynamic  and  complex  modulus 
approaches  is  characterized  and  the  equivalence  of  both 
approaches  is  proven  in  Chapter  4.  In  Chapter  5,  the 
damping  and  complex  stiffnesses  of  general   laminates  are 


derived  by  using  the  laminated  plate  theory,  the  energy 
approach,  and  the  preceding  derivations.  The  complex 
stiffnesses  are  completely  expressed  in  Appendix  A. 

The  environmental  conditioning  of  the  test  specimens, 
the  static  flexure  test,  and  the  impulse  hammer  techniques 
are  presented  in  Chapter  6.  These  experimental  methods, 
although  simple,  are  very  versatile  and  are  adequate  in 
determining  the  necessary  data  for  the  purpose  of 
this  investigation. 

The  theoretical  and  experimental  results  are  given  in 
Chapters  7-10.  The  moisture  and  thermal  expansions  of 
composites  are  quantified  in  Chapter  7.  The  current 
experimental  results  and  data  and  conclusions  of  previous 
investigators  are  used  in  Chapter  8  to  model  the  complex 
moduli  of  epoxy  as  functions  of  temperature  and  moisture 
content.  In  Chapters  9  and  10,  the  hygrothermal  effects  on 
the  stress  field  across  laminates  and  on  damping  of 
composites  are  investigated  with  the  help  of  the  results  in 
the  preceding  Chapters.  The  Finite  Element  Method  (F.E.M.) 
used  in  determining  the  stresses  is  summarized  in 
Appendix  B. 


CHAPTER  2 
DIFFUSION  OF  MOISTURE 


2.1  Introduction 


The  mechanism  of  moisture  absorption  and  desorption 
in  most  fiber  reinforced  composites  is  adequately  described 
by  Fick's  law  [4].  Fick  recognized  that  heat  transfer  by 
conduction  is  analogous  to  the  diffusion  process. 
Therefore,  he  adopted  a  mathematical  formulation  similar  to 
Fourier's  heat  equation  to  quantify  the  diffusion 
process  [24 ,  25] . 

2 . 2  Fickian  Diffusion 

The  Fourier  and  Fick's  equations,  describing  the 
one-dimensional  temperature  and  moisture  concentration,  are 
respectively  given  by 
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where 


is  the  density  of   the  material 


i  s  the 


specific  heat,  T  is  the  temperature,  t  and  x  are  the  time 

and   spatial   coordinates,   respectively,   K    is   the   thermal 

conductivity,   c  is   the  moisture  concentration,   and  D    is 
J  x 

the  moisture  diffusivity. 

The   moisture   diffusivity,   D  ,   and   the   thermal 

x 

diffusivity,  K  /(pC  ),  are  the  rate  of  change  of  the 
moisture  concentration  and  the  temperature,  respectively. 
In  general,  both  parameters  depend  on  temperature  and 
moisture  concentration.  But  experimental  data  show  that, 
for  most  composites,  moisture  diffusivity  does  not  depend 
strongly  on  moisture  concentration  [4].  Hence,  Eq  (2.2) 
becomes 
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(2.3) 


and  is  solved  independently  of  Eq  (2.1). 

The   three-dimensional   diffusion   in   an   anisotropic 
medium  is  obtained  by  generalization  of  Eq  (2.2)  as  follows 
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where  the  diffusivity  matrix  is 
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Expansion  of  Eq  (2.4)  results  in  an  equation  of  the  form 
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if  the  coefficients  D.  .'s  are  considered  to  be  constant. 

i  J 


2 . 3  Fickian  Diffusion  in  a  Plate 

Laminated  plates  are  widely  used  in  the  experimental 
characterization  of  composites.  Hence,  being  of  practical 
interest,  the  problem  of  moisture  absorption  in  a  plate  is 
thoroughly  discussed  in  this  section. 

2.3.1  Infinite  Plate 

The  case  of  moisture  absorption  through  a  material 
bounded  by  two  parallel  planes  is  considered.  The  initial 
and  boundary  conditions  of  an  infinite  plate  exposed  on 
both  sides  to  the  same  constant  environment  (Figure  2.1) 
are  given  by 
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T  =  T. 
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for  z  =  0,  z  =  h  and  t  >  0 


where   T.   is   a   constant   temperature,   c.   is   the   initial 
1  l 

moisture  concentration  inside  the  material,  and   c    is  the 

00 

maximum  moisture  concentration.  It  is  assumed  that  the 
moisture  concentration  on  the  exposed  sides  of  the  plate 
reaches   c    instantly. 

00  J 

The   solution   of   Eq  (2.3)   in   conjunction  with   the 
conditions  of  Eqs .  (2.7)  is  given  by  Jost  [25] 
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Equation  (2.8)  is  plotted  in  Figure  2.2. 

The  average  moisture  concentration  is  given  by 
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Substitution  of   Eq .  (2.8)   into  Eq .  (2.9)   and   integration 
result  in 
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This  analysis  can  be  applied  to  the  case  of  diffusion  of 
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moisture   into   a   laminated   composite   plate   so   thin   that 
moisture  enters  predominantly  through  the  plane  faces. 

2.3.2   Semi -Inf ini te  Medium 

In  the  early  stages  of  moisture  diffusion  into  a 
plate,  there  is  no  interaction  between  moisture  entering 
through  different  faces.  Therefore,  the  solution  of 
moisture  absorption  into  a  semi- inf ini te  half-plane  is 
applicable  to  a  plate  for  short  time. 

The  initial  and  boundary  conditions  of  a  semi- 
infinite  plane  (Figure  2.3)  exposed  to  a  constant  moist 
environment  are 
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for   z  =  0  and  t  >  o 


The  solution  of  Eq  (2.3)  in  this  case  is  [24,  25] 
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The  rate  at  which  the  total  specific  mass  of  moisture,  m, 
is  diffusing  into  the  half-plane  is 
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Thus,  the  total  mass  of  moisture  entering  through  an  area  A 
in  t  ime  t  is 
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Equation  (2.14)  shows  that  the  mass  of  diffusing  substance 
is  proportional  to  the  square  root  of  time. 

2.3.3   Experimental  Measurement  of  Moisture  Content 

In  the  case  of  a  finite  plate,   the  total  moisture 
con  tent  i  s 


m  =  pVc 


(2.15) 


where  V  is  the  volume  of  the  test  piece.  The  total  moisture 
content  is  experimentally  measured  by  subtracting  the  dry 
weight,  w  ,  ,  from  the  current  weight,  w,  of  the  plate,  i.e, 


m  =  w  -  w 


(2.16) 


A  parameter  of  practical  interest  is  the  percent  moisture 
content  defined  as 
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Since 
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the  experimentally  measured  M  of  Eq  (2.18)  can  be  compared 
to  the  analytical  value  given  by  Eq  (2.10). 

2.3.4   Approximate  Solution  of  Moisture  Content 

Approximate  solutions  of  the  specific  moisture 
distribution  in  a  plate  subjected  to  the  conditions  given 
by  Eqs  (2.7)  are  useful,  since  the  difficulty  of  dealing 
with  infinite  series  can  be  avoided. 

Smal 1  t  ime .  As  discussed  in  section  2.3.2,  Eq  (2.14) 
can  be  applied  during  the  early  stages  of  absorption.   It 
yields 
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Large  t  ime .  Tsai  and  Hahn  [26,  p.  338]  suggest  that, 
for  sufficiently  large  t,  Eq .  (2.10)  can  be  approximated  by 
using  the  first  term  of  the  series,  i.e., 


13 


C    Ci    ,    8 

=  1  -  — ^-exp 


c™~  c  • 


U 


-  IT 


2Dzt 


(2.20) 


Shen  and   Springer   formulation.     These   researchers 
have  derived  in  Ref.  [4]  the  following  approximation 
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Figure  2.4  shows  a  comparison  of  Eqs  (2.19-21)   with  the 
exac  t  so lu t  i  on . 

2.3.5   Edge-Effect  Corrections  in  a  Finite  Laminated  Plate 

A  plate  exposed  to  a  humid  environment  absorbs 
moisture  through  all  its  six  sides.  At  small  time,  the 
interaction  of  moisture  entering  through  different  sides  is 
negligible.  Therefore,  Eq .  (2.19)  can  be  applied  to  such 
cases .  It  yields 


m  =  4P(Coo"  ci) 
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where  D  ,  D  ,  and  D   are  the  di f f us i vi t ies  in  the  x,  y,  and 
x    y        z 

z  directions,  respectively.  The  geometry  of  the  plate  is 
shown  in  Figure  2.5.  Rewriting  Eq .  (2.22)  in  terms  of  the 
percent  moisture  content  gives 
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where  the  effective  diffusivity  D  is 
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The  mi cromechani cs  formulation  for  di f f us ivi t ies  proposed 
by  Shen  and  Springer  [4]  and  modified  by  Hahn  [26]  for 
impermeable,  circular  cross-section,  fiber-reinforced 
compos  i  tes  is 
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where   D  .   D_   and   D,   are   the   matrix,   transverse,   and 
m     T         L 

longitudinal  di f f us i vi t ies ,  respectively.  Equation  (2.23b) 
for  a  unidirectional  lamina  with  all  fibers  parallel  to  the 
x-direction  can  be  written  as 
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For  a  general  laminated  plate  consisting  of  N  layers  with 

fiber  orientations   6.,    the  di f f us i v i t i es  are 
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where  h.  is  the  thickness  of  the  j-th  layer.  The  effective 
J 

diffusivity   of   a   general    laminate    is    obtained   by 
substituting  Eqs.  (2.26)  into  Eq  (2.23b). 


2 . 4  Diffusivity  and  Maximum  Moisture  Content 


The  diffusivity  D  and  the  maximum  moisture  content 
M  must  be  experimentally  determined  in  order  to  predict 
the  moisture  content  and  distribution  in  a  lamina.  These 
parameters  are  obtained  by  the  following  procedures: 

-  a  thin,  unidirectional  composite  plate  is 
completely  dried  and  its  weight  is  recorded. 
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-  the  specimen  is  then  placed  in  a  constant 
temperature  and  constant  relative  humidity  environment,  and 
its  weight  as  function  of  time  is  recorded, 

-  the  moisture  content,  M,  versus  the  square  root  of 


time,    /t ,   is  plotted  as  shown  in  Figure  2.6. 

The  maximum  moisture  content  is  determined  from  the 
plot  and  the  diffusivity  from  the  following  equation 


w2 :  Mi 


=  4M. 
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The  subscripts  1  and  2  are  defined  in  Figure  2.6. 

The   diffusivity   depends   only   on   the   material   and 
temperature  as  follows 


D   =  D  exp 
z     o 


RT 


(2.28) 


where   R   is   the   gas   constant,   D    and   E.   are   the 

o  d 

pr e-exponen t ial    factor    and    the    activation    energy, 
respect  i ve ly . 

Experimental  research  has  shown  that  the  maximum 
moisture  content  depends  on  environment  humidity  content 
and  material.  For  a  material  exposed  to  humid  air  [4],  the 
equilibrium  moisture  content  can  be  expressed  as 


Mm  =  a<J>1 


(2.26) 
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where   (j)   is  the  relative  humidity,  a  aad  b  are  material 
cons  tan t  s . 


IS 


Moisture 


Moisture 


x* 


Fig.  2.1   Plate  subjected  to  a  constant  humid  environment 
on  both  sides. 
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Fig.  2.2    Moisture   distribution   across   a   plate.   The 
numbers   on   the   curves   are   the   values   of 
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Fig.  2.3    Semi-infinite  plate  in  a  humid  environment 


—  Exact 
-  One— term 
•    Shen  and  Springer 


Fig.  2.4  Comparison  of  the  exact  specific  moisture 
concentration  equation  with  some  approximate 
so lu  t  ions . 
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Fiber  direction 


Fig.  2.5    Geometry  of  a  plate 
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Fig.  2.6    Moisture  content  versus  square  root  of  time.  On 
the  curve   /  t   <   /  t_  <   /  t.  and  the  slope  is 

constant  for  /t  <   /  t.  . 


CHAPTER  3 
COMPLEX  MODULI  OF  UNIDIRECTIONAL  COMPOSITES 


3 . 1  In  t roduc  t  i  on 

Composite  materials,  such  as  Glass/Epoxy  and 
Graphi te/Epoxy ,  have  a  polymeric  matrix.  Therefore,  they 
display  viscoelastic  behavior.  Some  of  the  effects  of  this 
time-dependent  phenomenon  are:  stress  relaxation  under 
constant  deformation,  creep  under  constant  load,  damping  of 
dynamic  response,  etc. 

This  chapter  is  an  introduction  to  the  dynamic 
behavior  of  viscoelastic  composites  in  terms  of  complex 
modu  1  i  . 

3 . 2  General  Theory 

A  usual  representation  of  the  one-dimensional 
stress-strain  relation  of  a  viscoelastic  material  subjected 
to  a  harmonic  strain  history  of  the  form 


€(t)  =  €°eiCOt  (3.1) 
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is  g  i ven  by 


a(t)  =  B*(iw)€°e1(Jt  =  B*(iu)€(t) 


(3.2) 


The  complex  modulus  B   can  be  decomposed  into  its  real  and 
imaginary  parts  as  follows 


B  (iu)  =  B'  (w)  +  iB"((j) 


(3.3) 


The  terms  B'  and  B"  are  called  the  storage  and  loss 
moduli,  respectively,  and  the  ratio  of  the  loss  over  the 
storage  modulus 


n  = 


IT 

B' 


(3.4) 


is  referred  to  as  either  the  loss  factor  or  damping.  The 
loss  modulus  is  a  measure  of  the  energy  dissipated  or  lost 
as  heat  per  cycle  of  harmonic  deformation. 

3 . 3  Mi cromechani cs  Formulation  of  Elastic  Moduli 


The  longitudinal  modulus  E ..  1  ,  the  transverse  modulus 
E„„,  the  in-plane  shear  modulus  G  „,  and  the  major 
Poisson's  ratio  v  „  can  be  obtained  by  using  the  rule  of 
mixtures  and  the  Halpin-Tsai  equations,  viz., 


E,   =  v.E     +  v  E 
11     ffll     mm 


(3.5) 
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1  +  2n  v. 

p     —  F   
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G,„  =  G 


1  +  n2vf 


12     ml-  n„v. 


(3.7) 


12=   ffl2     mm 
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where 


n,  = 


(Ef22/Em>  -  1 
(Ef22/Em)  +  2 


(3.9) 


(Gfl2/Gm)  -  1 
"2  =  ^Gfl2/Gm)  +  l 


(3.10) 


The   subscripts   f   and   m   stand   for   fiber   and   matrix 

respectively,  and  vr  and  v   are  the  volume  fractions. 

r        m 


3 . 4  Comp 1  ex  Modu 1 i 


The  mi cromechani cs  formulations  of  the  complex  moduli 
are  obtained  by  applying  the  e las t ic-vi scoe las t i c 
correspondence  principle  [27-29],  i.e.  ,  by  undertaking  the 
following  steps: 

i)  determining  the  elastic  moduli  of  composites  in 
terms  of  the  constituent  material  properties, 
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ii)  replacing  the  elastic  ~  o  cl  u  1  i  o :  fibers  a  r.  cl  -  a  t  r  i  x 
by  corresponding  complex  expressions. 

For   a   viscoelastic   composite,    the  properties   o :    the 
constituent  materials  are 


Ef 1 1  "  Ef 1 1  +  lEfll 
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G   = 


G'  +  iG" 
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v      = 


V        +  1  u 
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Dfl2  =  Dfl2 


The   bulk   modulus   of   eooxv   matrix.   K  .    is   real   and 

m 

independent  of  frequency  [2S].  It  is  given  by 


m  "  3(1  -  2d  ) 


(3.12) 


while  the  viscoelastic  bulk  -.ecu  1  us   is  obtained  from  the 
correspondence  principle 
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Separation     of      the     real     and     imaginary     parts     of     Eq .     (3.13) 
yields 
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(3.14) 
Since  the  dilatation  bulk  modulus  is  real,  the  imaginary 
part  of  Eq .  (3.14)  is  equal  to  zero;  hence 


2E"l>"  +  E"(l  -  2v  •  )  =  0 
mm     m v        m' 


(3.15) 


Equation  (3.15)  results  in 
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(3.16) 


The  shear  modulus  of  the  matrix  is  given  by 


G*  = 
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m 


2( 1  +  v    ) 


3K    E^ 
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9K 
m 
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m 
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Separating  the  real  and  imaginary  parts  and  neglecting  the 

2 
terms  of  the  form   (E")   yield 

v  m'  J 
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G       = 
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9K  E" 


1    +    i 


m  m 


9K       -    E"     E' 
m  mm 


(3. 17b) 


Introduction    of     the    material    properties 
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into    Eqs.     (3.11)     results     in 


Efn     =    Ef  1  A1    +    iT1<-n  ) 


f  1 1    ~       f  1  1 


f  11 


fOO     ~         f  OOV  +      1T)flO' 


f22    "       f22 


f22 


f  12 


Gfl2(1    +    iT?fl2> 


E       =    E" (1    +    in    ) 
m  m v  m 


(3.19) 


G      s   G'{1    +    ill.    ) 
m  mv  'Gra; 
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v       = 
m 


v  '  +  in  ( v ' 
m      mm 


f  12 


Ufl2  "  Ufl2 


There  are  no  satisfactory  data  on  the  shear  and  transverse 

damping  of   fibers.   Fibers  have  damping  with  a  magnitude 

order   ten   times   smaller   than  epoxy .   The  dampings    n    , 

Tjf<?9,  and  n .  ^   are  assumed  to  be  equal  and  are  replaced  by 

n    in  subsequent  equations.  Since  the  fiber  damping,   n 

is   much   smaller   than   the   matrix   damping,     n  ,     the 

m 

imaginary  part  of  the  fiber  Poisson's  ratio  is  neglected. 
The  preceding  assumptions  have  a  negligible  effect  on  the 
complex  moduli  of  composites. 

Application  of  the  e las t i c-vi scoe las t i c 
correspondence  principle  to  Eqs.  (3.5)-(3.10)  and 
substituting  them  into  Eqs  (3.19)  yield  the  following 
complex  material  properties 


E    =  v.E'   (1  +  ii7  )  +  v  E'(l  +  in  ) 
11     ffllv        f '  mmv        m' 


1  +  2n  v 

E99  =  Em(l  +  in) -^ 

2.2.  m        m  ,     »* 

1  -  nlVf 


(3.20) 


x  l    +    n2Vf 

G-   :  Gm(l  +  inGm)- !-£■ 

1  -  n2vf 
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where 


12 


(v.i).,„  +  v  v  '  )  +  i  T7  (  u  ' 
f  f  1  2     mm       mm 


bvm 


*         E'22(l  +  inf)  -  Em(l  +  iVJ 

n   =  

E-  0(1  +  iT7f)  +  2E'(1  +  iT7m) 


f22 


m 


m 


n2  = 


G'12(l  +  iT?f)  -  Gm(l  +  ir,^) 
Gfl2(l  +  It,,)  +  G;(l  +  i„r_J 


The  elastic  moduli  given  by  Eqs  .  (3.5)-(3.8)  model 
experimental  results  with  a  good  accuracy  [2].  Therefore, 
they  are  used  instead  of  mathematically  exact 
micromechani cs  formulas,  such  as  those  derived  by 
Hashin  [27,  29]. 


CHAPTER  4 
DAMPING 


4 . 1  Damping  Mechanisms 

Any  vibrational  energy  introduced  in  a  structure 
tends  to  decay  in  time.  This  phenomenon  is  called  damping. 
There  are  two  types  of  damping  mechanisms,  external  or 
nonmaterial  and  internal  or  material. 

4.1.1  Nonmaterial  Damping. 

Two  common  types  of  external  damping  are 

-  Accoustic  damping:  a  vibrating  structure  always 
interacts  with  the  surrounding  fluid  medium  (air,  water, 
etc.).  This  effect  can  lead  to  noise  emission  and  even  to 
changes  of  the  natural  frequencies  and  mode  shapes.  Thus, 
mechanical  responses  might  be  modified. 

-  Coulomb  friction  damping:  two  contacting  surfaces 
in  relative  motion  dissipate  energy  through  frictional 
forces . 
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4.1.2  Material  Damping 

There  are  many  damping  mechanisms  that  dissipate 
vibrational  energy  inside  the  volume  of  a  material.  Damping 
phenomena  include  thermal  effects,  magnetic  effects,  stress 
relaxation,  phase  processes  in  solid  solutions  [30,  p.  61], 
etc. 

The  internal  damping  of  polymeric  matrix  composites, 
such  as  Glass/Epoxy  and  Graphi t e/Epoxy ,  is  dominated  by 
viscoelastic  damping. 

4 . 2  Characterization  of  Damping 

4.2.1  Free  Vibrati  on 

A  cantilever  under  free  vibration  oscillates 
regularly  with  an  amplitude  that  decreases  from  one 
oscillation  to  the  next  one  (Figure  4.1).  A  measure  of 
damping  is  the  logarithmic  decrement  defined  as 


r  A 


6  =  Nln 


uAn+NJ 


(4.1) 


where 

A   =  amplitude  of  the  n-th  cycle 

n      K  J 

A   AI  =  amplitude  of  the  (n+N)-th  cycle 
n+N      r  y  j  j 

The  damping  defined  in  Eq  (4.1)   is  applicable  to  viscous 


31 

damping  and  for  hysteretic  damping  that  is  represented  by  a 
complex  modulus  approach. 

4.2.2  Steady  State  Vibration 

Damping  also  influences  the  dynamic  equilibrium 
amplitude  of  structures  (e.g.  beams)  that  undergo  harmonic 
oscillation.  A  resonance  usually  occurs  (Figure  4.2).  The 
following  measure  of  damping  is  used 


o   -  w 

V    =  — L  (4.2) 

o 


where 


cj   =   resonant  frequency 
o  M 


w1  ,  w~  =   frequencies  on  either  sides  of  w    such  that 


the  amplitude  is    1/  /2     times  the  resonant 
amp  1 i  tude . 

In  the  case  of  a  vibration  induced  by  the  force 

f ( t )  =  Fsin(ut) 

the  response  (deflection),  w(t),  is  out  of  phase  with  f(t) 
by  an  angle  e.       such  that 

w( t )  =  Wsin(wt  +  e) 
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The  work  done  per  cycle  is 


n2I7/cj 


D  = 


i(t)p-   dt 
v  'dt 


=  HWF  sin(e) 


(4.3) 


^o 


The  strain  energy  stored  in  the  system  at  the  maximum 
displacement  is  half  the  product  of  the  maximum 
displacement  by  the  corresponding  value  of  the  force,  i.e.. 


U  =  ^FW  cos(fe) 


(4.4) 


There  is  no  damping  if  the  work  done  per  cycle  is  zero, 
i.e.  if  sin(e)  =  0. 

The  ratio  of  energy  dissipated  in  a  cycle  to  energy 
stored  at  the  maximum  displacement  is  another  measure  of 
damping.  Therefore,  the  damping  is 


77  =  2lU  =  tan{£) 


(4.5) 


The  definitions  of  damping  given  by  Eqs.  (4.2)  and  (4.5) 
are  equivalent  [33]. 

4.2.3  Complex  Modulus  Approach 


The  one-dimensional  stress-strain  relation  of  a 
viscoelastic  material  undergoing  harmonic  motion  has  been 
shown  to  be  (Eq.  (3.2)) 
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a*(t)  =  E*(o)€o  e1Cl,t  =  (E'(u)  +  iE"(u))€   eiwt   (4.6) 


Noting  that   i  |  (J  |  €  =  d€/dt,   Eq .  (4.6)  can  be  written  as 


*r  +  \  c       lUt  ,   E   .   iut 

a    (t)  =  €  e     +  -i — r    lwe    € 
v  '  O  GJ  o 


(4.7) 


The  real  part  is  given  by  (after  algebraic  manipulation) 


o( t )  =  E '€  si 

v  '  o 


in(wt  +  e)  -i  1  +  V 


(4.8) 


where   n  =  tan(e.)  =  E"/E' 

The  energy  dissipated  during  a  cycle  per  unit  volume  is 


D  =  ()  a    d€   = 


■217/u 


d€ 

°    ~dt  dt  = 


HnE'62 
o 


(4.9) 


The  maximum  energy  stored  is 


1  2 

U    =    k   E'     €« 

2  o 


(4.10) 


Therefore , 


r\    = 


LI         _2_ 
E"     ::    2/TU 


(4.11) 
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Hence,  the  definitions  of  damping  given  by  Eq .  (3.4)  and 
Eq .  (4.5)  are  equivalent.  This  conclusion  is  also  valid  for 
general  cases  of  structural  vibration. 
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Envelope  Equation 


Vibrating  beam 


time 


Response 


Fig.  4.1  Schematic  drawing  of  a  free-clamped  beam  under 
free  vibration  and  plot  of  its  deflection 
versus  t  ime . 


Fe 


lot, 


Beam  under  forced  vibration 

•A 


No  damping 


Deflection  versus  time 


Amplitude  versus  frequency 


Fig.  4.2  Schematic  drawing  of  a  free-clamped  beam  under 
forced  vibration  and  plots  of  the  deflection 
versus  time  and  deflection  amplitude  versus 
frequency . 


CHAPTER  5 
DAMPING  AND  STIFFNESSES  OF  GENERAL  LAMINATES 


5  .  1  In t  roduc  t  i  on 

Both  the  laminated  plate  theory  and  the  energy  method 
approaches  for  analyzing  the  damping  and  the  stiffnesses  of 
general  laminates  are  presented  in  this  chapter. 

5 . 2  Laminate  Plate  Theory  Approach 

Four  independent  parameters  are  needed  to  determine 

completely  the  damping  of  a  unidirectional  composite.  But, 

the  analysis  of  the  material  damping  of  a  general  laminated 

composite  requires   the  use  of   eighteen  parameters.   These 

quantities  are  the  ratios  of  the  imaginary  over  the  real 

parts   of   the   complex   in-plane   stiffnesses   A. .'s,   the 

complex  coupling  stiffnesses  B.  .'s,  and  the  complex  bending 

stiffnesses  D. .'s  (see  Appendix  A), 
ij 

The  terms  A.  .'s,  B.  .'s,  and  D.  .'s  are  defined  as 


A*. 


■+h/2 


-h/2 


Q.  .  dz 
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B*.  = 


■+h/2 


-h/2 


zQ.  .  dz 


(5.1) 


D 


ij 


r+h/2 


-h/2 


z  Q.  .  dz 


where  the  complex  transformed  stiffness    Q.  .'s    depend  on 

ij 
x     m     ^      m 
E .,  ,  .  E„~,  G  „,    v-\o       and  the  orientation  of  each  layer  of 

the  laminate . 

The  in-plane,  coupling,  and  flexural  material  damping 

are  defined  as 


InU 


i  J 


C^ij 


BV  . 
i  J 

b:  . 

ij 


(5.2) 


DV  . 
i  J 

d:  . 


respec  t  i ve ly 


5 . 3  Energy  Method  Approach 


The  energy  method  can  be  used  to  determine  the 
damping  of  laminated  composite  materials  under  certain 
loading  and  boundary  conditions.  The  damping  of  a  laminated 


3S 


composite  material   in  the  first  mode  of  vibration  can  be 
de  f  ined  as 


N 


1    (kUd> 


eye 


I  = 


k=l 


(5.3) 


2U  .  U 
k  s 


k=l 


where  N  is   the   total  number  of   layers,   (.  U.)       is   the 

J  vk  d' eye  . 

energy  dissipated  in  the  k-th  layer  during  a  cycle,  and  ,U 
is  the  maximum  energy  stored  in  the  k-th  layer.  The  storage 
and  the  dissipated  energy  are  given  by 


kUs  =  I 


e.c: .€.  dV 

J  Jl  1 


U  ,  =  IT     £  .C"  .€. 

'-    d      Jv   j  ji  l 


(5.4) 


dV 


where  i  and  j  are  the  material  principal  axes,    C '.  .    and 

J  i 

C"    are   the   real   and   imaginary   parts   of   the   complex 

J  i 

stiffnesses  and  V,   is  the  volume  of  each  layer.  Hence,  the 

k  J 

"total"  damping  of  an  N-layered  laminate  is  given  by 


N 


}    J   {€}T[C"]{€}  dV 


n  = 


k=l  Vk 


l    J   {€}T[C']{€}  dV 


(5.5) 


k=l   k 
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The  maximum  strain  vector  {€}  can  be  determined  by  the 
finite  element  method  first.  Then,  the  damping  can  be 
deduced.  Equation  (5.5)  is  used  to  determine  the  damping  of 
a  beam  with  variable  thickness  or  of  more  general 
s  true  tures . 


CHAPTER  6 
EXPERIMENTAL  PROCEDURES 


6  .  1  In  t  r oduc  t  i  on 


A  description  of  the  test  specimens  and  the 
experimental  procedures  of  the  present  investigation  is 
given  in  this  chapter. 

6 . 2  Test  Specimens 

The  test  specimens  used  to  determine  the  complex 
moduli  of  epoxy  and  of  composite  materials  are  thin  strips 
of  approximate  dimensions,  150mm  by  25mm  by  2mm.  The  only 
materials  tested  are  Magnolia  2026  laminating  epoxy  and  3M 
Scotchply  Glass/Epoxy.  The  curing  temperatures  of  the  epoxy 
and  the  Glass/Epoxy  are  175  C  and  170  C,  respectively.  The 
initial  properties  of  these  materials  (at  20  C  and  without 
moisture),  as  well  as  those  of  a  typical  Graphi te/Epoxy , 
are  given  in  Table  6.1. 
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6 . 3  Environmental  Conditioning 

The  specimens  are  conditioned  in  a  Thermotron 
environment  chamber  at  a  constant  temperature  and  constant 
relative  humidity.  The  weight  gain  of  the  test  pieces  as  a 
function  of  time  is  monitored.  Right  after  moisture 
equilibrium  is  reached,  the  specimens  undergo  all  tests  at 
diverse  temperatures  inside  a  testing  chamber  connected  to 
the  environment  chamber  (Figure  6.1).  The  range  of 
temperature  achieved  inside  the  environment  chamber  is  4  C 
to  90  C  and  the  range  of  relative  humidity  is  4%  to  99%  for 
temperatures  below  75  C.  As  temperature  increases,  the 
highest  relative  humidity  that  can  be  obtained  decreases 
steadily  to  75%  at  90°C. 

6 . 4  Four-Point  Flexure  Test  Method 

The  Young's  modulus  and  the  Poisson's  ratio  can  be 
determined  with  the  four-point  flexure  test  method.  The 
loading  configuration  of  this  test  is  shown  in  Figure  6.2. 
The  elastic  flexural  analysis  yields  [31] 


E  = 


PI 


8bh3w 


(6.1) 


where   E   is  the  effective  modulus,  P  is  the  applied  load, 
1  is  the  length  of  the  specimen,  b  is  the  specimen  width,  h 
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is  the  thickness,  and  w  is  the  deflection  at  quarter-point. 
Poisson's  ratio  is  expressed  as 


v    = 


6hw   y 


(6.2) 


where   the   transverse   strain    €      is   measured   with   a 

y 

transverse   strain   gage   cemented   in   the   middle   of   the 
spec  imen . 

6 . 5  Impulse  Hammer  Technique 


The  material  damping  and  the  storage  modulus  of  a 
one-dimensional  thin  beam  are  determined  with  the  impulse 
hammer  technique.  This  technique  was  pioneered  by  Halvorsen 
and  Brown  [32].  The  equipment  set-up  is  shown  in 
Figure  6.3.  The  specimen  is  clamped  inside  the  testing 
chamber.  A  force  impulse  is  applied  to  the  test  piece  by  a 
force  transducer.  The  end  displacement  of  the  specimen  is 
recorded  with  a  non-contacting  motion  transducer.  Both 
responses  from  the  force  and  motion  transducers  go  through 
signal  conditioning  equipments  (filters,  amplifiers).  These 
responses  are  digitized  in  a  Fast  Fourier  Transform 
analyzer  (FFT)  to  obtain  the  transfer  function  in  terms  of 
the  frequency.  The  transfer  function  is  defined  as  the 
ratio  of  the  Fourier  Transform  of  the  output  (displacement 
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v(t))   over   the   Fourier   Transform   of   the   input   (force 
impulse  u(t));  that  is, 


H(f)  -  ^ill 
nU  ;    U(f  ) 


(6.3) 


where 

t  =  t  ime 
f  =  frequency 
V(f)  =  Fourier  Transform  of  v(t) 
U(f)  =  Fourier  Transform  of  u(t) 
The  real  and  imaginary  parts  of  H(f)  are  displayed  on  the 
FFT  analyser  CRT  (Figure  6.3).  The  material  damping  defined 
by  Eq .  (4.11)  is  experimentally  obtained  by  the  following 
express  ion 


T7  = 


(f  /f,  )   -  1 

v  a  b' 

(f  /fK)2  +  1 
v  a  b' 


(6.4) 


where  the  frequencies  f   and  f,  are  defined  in  Figure  6.3. 
n  a       b  & 

The  storage  modulus  is  expressed  as  [33.  p. 464] 


E'  =  38.32  f2  p^r 


(6.5) 


where  f  is  the  resonant  frequency  in  Hz.  ,  p,  is  the 
material  density,  1  is  the  length  of  the  specimen  and  h  is 
the  thickness  of  the  specimen.  Equation  (6.5)  is  valid  only 
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for  the  case  of  the  first  mode  free  vibration  of  a 
c lamped- f ree  beam.  A  complete  description  and  analysis  of 
the  impulse  hammer  technique  are  presented  in  Lee's 
dissertation  [34]. 
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Table  6.1  Initial  properties  of  Magnolia  2026  epoxy, 
3M  Scotchply  Glass/Epoxy,  and  a  typical 
Graphi te/Epoxy  composite. 


Proper  ties 


Epoxy 


Glass/Epoxy     Graphi te/Epoxy 


0.50 


0.70 


p    (g/cnT 

) 

1  .25 

1  .93 

1.6 

En     (GPa. 

) 

4.0 

37.00 

155.23 

E22    (GPa. 

) 

4.0 

11  .54 

10.81 

G12    (GPa. 

) 

1.52 

3.46 

4.35 

U12 

0.32 

0.285 

0.217 

^11 

0.018 

0.0023 

0.0019 

^22 

0.018 

0.015 

0.0078 
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Testing  chamber 


Environmental  chamber 
/ 


MTS  load 
frame 


Fig.  6.1    Schematic  drawing  of  environmental  and  testing 
chambers . 


1/4 


PA 


1/4 


P/2 


P/2 


Fig.  6.2    Loading     configuration     of 
f 1 exure  test. 


the     4-point 
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Teat  specimen 
Hammed 


Clamp 


Motion  transducer 
•Testing  chamber 


Signal 

conditioning  equipment 


CRT 
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Fourier  Transform 


Real  part 


Im.  part 


fr      Frequency 


Fig.  6.3  Schematic  drawing  of  the  impulse  hammer 
technique  apparatus  and  a  typical  display  of 
the    Transfer    Fourier    Transform. 


CHAPTER  7 
HYGROTHERMAL  EXPANSION 


7.1  Introduction 


When  a  metallic  or  composite  structure  is  subjected 
to  a  change  of  temperature,  there  are  dimensional 
variations  and  there  may  be  stress  development.  For  a 
one-dimensional  case,  it  is  assumed  that  the  thermal  strain 
is  given  by 


€T  =  a. (T  -  T  )  =  a. AT  (7.1) 

1      l  v       o '      1  v     ' 


where 


a.  =  coefficient  of  thermal  expansion 
l 


T  =  actual  temperature 

T   =  reference  temperature. 
o  K 

A  polymer  matrix  composite  exposed  to  a  humid 
environment  absorbs  moisture.  Hence,  it  increases  in  weight 
and  dimensions.  This  situation  produces  a  moisture  strain 
that  varies  linearly  with  moisture  concentration  [26].  In 
the  one-dimensional  case,  the  hygrostrain  is  given  by 
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€H  = 


Pjfc 


c  ) 

o ' 


0.  Ac 
1 


(7.2) 


where  c   is  the  initial  moisture  concentration  and  0.    is 
o  1 

the  coefficient  of  moisture  expansion. 


7 . 2  Coefficients  of  Thermal  Expansion 

In   the   case   of   laminated   composite  plates,   three 

coefficients  of  thermal  expansion  are  used  in  determining 

the   thermal   strains.   These  parameters   can  be  written   in 

terms  of  fiber  and  matrix  properties.   The  mi cromechanics 
formulas  for  a  unidirectional  orthotropic  lamina  are  given 

by  (see  Refs.  [35.  p.  24],  and  [36,  p.  405]  for  a  detailed 
der  i va t  i  on) 


al  = 


v„arE.  +  v  a  E 

f  f  f m  m  m 

Ell 


<xn    =    v.a.    +    v    a       +    v.u.a.    +    v    v    a       -    u1_a1  (7.3) 

2  ff  mm  fff  mmm  12    1  v  ' 


«12  =  o 


where  the  subscripts  1  and  2  represent  the  fiber  and  the 
transverse  directions.  The  thermal  expansion  coefficients 
of  an  orthotropic  lamina  whose  fibers  make  an  angle  0 
with  the  x-direction  (Figure  2.5)  are  given  by 
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2  2 

a       =    a,     cos   0  +  a„  sin   0 
x     1  2 


2  2 

a   =  a,  sin   6  +  a_  cos   0 
y     1  2 


(7.4) 


a    =  2(a,  -  a~)  cos  0  sin  0 
xy     v  1     2y 


7 . 3  Coefficients  of  Moisture  Expansion 

Similarly,  the  coefficients  of  moisture  expansions  of 
an  orthotropic  lamina  with  impermeable  fibers  can  be 
expressed  as  [36,  p.  406] 


sE 


Pi     = 


m 


-P. 


1  s    E ,  ,     m 

mil 


2  sv  rn'm  121 


(7.5) 


m 


P12    =   o 


where  s  and  s   are  the  specific  gravities  of  the  composite 
m 

material  and  the  matrix.  The  moisture  expansion 
coefficients  expressed  in  an  axis  system  such  that  the 
x-direction  makes  an  angle  0  with  the  fibers  are  given  by 
Eqs.  (7.4)  after  replacement  of  a.'s    by  j3 .  's. 
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7 . 4  Experimental  Data 

7.4.1  Previous  Investigations 

Hahn  and  coworkers'  investigations  [21,  26,  37]  of 
swelling  of  composites  are  outlined  in  this  section.  Some 
of  the  typical  results  of  the  transverse  strain  versus 
percent  moisture  gain  are  obtained  by  conducting  the 
following  tests:  absorption  is  conducted  in  moisture 
saturated  air  such  that  Eqs  (2.2)  and  (2.7)  are  satisfied; 
while  desorption  takes  place  in  vacuum  at  the  same 
temperature.  Their  data  show  a  hysteretic  nature  of 
swelling  in  this  case. 

But,  when  swelling  of  composites  is  given  in  terms  of 
moisture  concentration,  the  average  behavior  of 
S2-Glass/Epoxy ,  Kevlar  49/Epoxy  and  Graphi te/Epoxy  can  be 
approximated  by 


e^  *  0.43c  =  P2c  (7.6) 


Since  the  data  presented  in  their  publications  display  a 
wide  scatter,  Hahn  et  al.  suggest  that  Eq .  (7.6)  can  be 
used  to  estimate  the  moisture  strains  for  most  composite 
mater  ial s . 


52 
7.4.2  Present  Investigation 

Epoxy  and  Glass/Epoxy  specimens  are  conditioned  at  a 
constant  relative  humidity  until  the  absorbed  moisture 
reaches  equilibrium.  The  changes  in  transverse  dimensions 
are  measured.  This  procedure  is  repeated  at  diverse  values 
of  relative  humidity.  The  results  are  plotted  in 
Figure  7.1.  The  longitudinal  swelling  strains  could  not  be 
measured  since  the  micrometer  calipers  used  were  not 
sufficiently  accurate.  These  data  yield  the  following 
experimental  values 


(3  (epoxy)  =  0.25 
m v       ' 


(7.7) 


02(Glass/Epoxy)  =  0.48 


(7.8) 


Substitution  of  Eq.  (7.7)  and  the  parameters  given  in 
Table  6.1  into  Eqs .  (7.5)  yields  the  following  empirical 
values 


P   =  0.042 


(7.9a) 


P2  =  0.47 


(7.9b) 


for   the   3M  Glass/Epoxy   composite.   The   experimental   and 
empirical  values  of   /3„  are  practically  equal.  Hence,   the 
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present  results  differ  slightly  from  the  approximation 
given  by  Eq .  (7.6). 

The  above  coefficients  and  the  typical  coefficients 
of  expansion  of  graphite  are  quantified  in  Table  7.1,  while 
the  storage  moduli  and  the  density  of  glass  and  graphite 
fibers  are  listed  in  Table  7.2.  These  properties  are  used 
to  plot  the  thermal  and  moisture  expansion  coefficients 
versus  the  fiber  volume  fraction  of  Glass/Epoxy  and 
Graphi te/Epoxy      in      Figures  7.2      through      7.5. 

The  values  in  these  plots  are  valid  for  dry 
composites  at  room  temperature.  Since  the  storage  modulus 
of  epoxy  varies  with  temperature  and  moisture  content,  this 
additional  effect  is  investigated  in  Chapter  8. 

In  general,  the  thermal  expansion  coefficients  are 
functions  of  temperature,  but  this  temperature  effect  is 
negligible  below  100  C.  Therefore,  in  the  subsequent 
chapters,  the  thermal  expansion  coefficients  are  assumed  to 
be  independent  of  temperature. 
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Table  7.1      Coefficients  of  moisture  and  thermal  expansion 
of  epoxy  and  graphite  and  glass  fibers. 


Epoxy 


Glass 


Graphi  te 


a    (um/m)/K     54.0 
P  0.25 


5.0 
0.0 


0.2 
0.0 


Table  7.2     Properties  of  Glass  and  graphite  Fibers. 


Glass 


Graphi  te 


E'n(GPa) 


Ef22(GPa) 


G^12(GPa) 


Vt 


f  12 


70.0 
70.0 
28.7 

0.0015 

0.22 


220.0 
16.6 
8.27 
0.0015 
0.  16 


p  (g/cm  ) 


2.60 


1  .75 
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Fig.     7.1         Transverse      moisture      strain     of      Magnolia     epoxy 
and    3M    Scotchply    Glass/epoxy. 
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Fig.  7.2    Plot  of  the  thermal  expansion  coefficients  in 
terms   of   fiber   volume   fraction   of   a   dry 

S  Glassf iber/Epoxy  at  20°C. 
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Fig.  7.3    Plot  of  the  thermal  expansion  coefficients  in 
terms   of   fiber   volume   fraction   of   a   dry 

Graphi te/Epoxy  at  20°C . 
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Fig.  7.4     Plot  of  the  moisture  expansion  coefficients  in 
terms   of   fiber   volume   fraction   of   a   dry 

S  Glassf iber/Epoxy  at  20°C. 
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Plot  of  the  moisture  expansion  coefficients  in 
terms   of   fiber   volume   fraction   of   a   dry 

Graphi te/Epoxy  at  20°C. 


CHAPTER  8 
HYGROTHERMAL  EFFECTS  ON  COMPOSITE  COMPLEX  MODULI 


8  .  1  Literature  Survey 

The  storage  moduli  (real  parts  of  Eqs.  (3.11))  of 
composites  are  usually  determined  by  dynamic  testings,  such 
as  the  technique  described  in  section  6.5.  They  can  be 
approximated  by  using  static  tests  [38]. 

Shen  and  Springer  [12]  investigated  the  environmental 
effects  on  the  elastic  moduli  of  a  Graphi te/Epoxy  composite 
and  made  a  survey  of  existing  data  showing  the  effects  of 
temperature  and  moisture  on  the  elastic  modulus  of  several 
composites.  Their  conclusions  are  listed  below. 

i)   The   hygrothermal   effects   on   the   0    fiber 
direction  laminates  are  negligible. 

ii)  For  90  fiber  direction  laminates,  the 
hygrothermal  effects  on  the  modulus  are  insignificant  in 
the  200K  to  300K  temperature  range.  But,  between  300K  and 
450K,  the  hygrothermal  effects  on  the  modulus  are 
irnpor  tan  t  . 

Putter  et  al.  [38]  investigated  the  influence  of 
frequency   and   environmental   conditions   on   the   dynamic 
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behavior   of   Graphi t e/Epoxy   composites.    Their   overall 
conclusions  are 

i)  The   effects   of   frequency   on   the   modulus   and 
damping  are  quite  small  in  all  cases. 

ii)  The  effects  of  frequency  on  the  modulus  and 
damping  are  relatively  greater  for  matrix-controlled 
laminates  at  higher  frequencies  (above  400  Hz.). 

iii)  At  the  same  temperature,  damping  increases  with 
moisture  saturation.  But  for  dry  laminates,  damping 
decreases  slightly  as  temperature  increases. 

From  all  these  experimental  works,  a  general  summary 
can  be  drawn:  the  influence  of  hygrothermal  conditions  on 
the  elastic  modulus,  dynamic  modulus  and  damping  of 
composites  is  matrix  dominated. 

8 . 2  Theoretical  and  Experimental  Assumptions 


Since  the  hygrothermal  influence  on  composite 
properties  is  matrix  controlled  [12,  38],  the  fiber 
properties  are  assumed  to  be  constant  at  any  temperature 
below  the  glass  transition  temperature  and  at  any  moisture 
content.  Therefore,  to  obtain  the  values  of  the  complex 
moduli  of  composites,  it  is  sufficient  to  know  how 
temperature  and  moisture  affect  the  complex  moduli  of  the 
epoxy  matrix,  and  then  use  the  mi cromechan i cs  formulations 
given  by  Eqs.  (3.20).  Thus,  only  the  following  functions 
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E"  =  E" (T.c) 
m     m 


v  '     -    v  '  ( T  ,  c  ) 
m     m 


(8.1) 


T?m  =  1  (T.c) 
m     m 


need  to  be  experimentally  evaluated.  The  constant  fiber 
properties  are  given  in  Table  7.2.  The  effects  of  frequency 
are  negligible  below  400  Hz.  The  results  of  this 
investigation  are  not  accurate  for  higher  frequencies  since 
their  effects  have  not  been  taken  into  account. 

The  qualitative  influence  of  temperature  only  on  the 
storage  modulus,  real  part  of  Poisson's  ratio,  and  damping 
of  epoxy  is  illustrated  in  Figures  8.1-8.3.  There  are  three 
distinct  regions.  At  room  temperature  (in  the  glassy 
region),  the  storage  modulus,  Poisson's  ratio,  and  damping 
of  epoxy  are  equal  to  about  4.0  GPa,  0.35,  and  0.018, 
respectively.  In  the  glassy  region,  the  storage  modulus 
decreases  slowly,  while  Poisson's  ratio  and  the  damping 
increase  as  temperature  increases.  In  the  next  region 
(transition  region),  the  storage  modulus  decreases  rapidly, 
and  both  Poisson's  ratio  and  damping  reach  their  maximum 
values.  The  last  region  is  the  rubbery  region  where  the 
modulus  takes  a  very  low  value,  and  all  three  parameters 
stay  relatively  constant. 

Typical  values  of  the  modulus  in  the  rubbery  region 

-2 

could  be  10   times  the  glassy  modulus  or  lower.  The  damping 
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can  reach  a  value  of  1  or  even  2  in  the  transition  region 
[30,  p.  90].  Poisson's  ratio  reaches  the  limiting  value  of 
0.5,  which  is  approximated  by  incompressible  rubbers 
[39,  p.  293]. 

The  position  of  the  transition  region  depends  on  the 
moisture  concentration.  The  effects  of  moisture  on  the 
glass  transition  temperature,  T  ,  of  six  epoxy  resins  have 
been  determined  by  Delasi  and  Whiteside  [6].  These  results 
are  plotted  in  Figure  8.4.  They  are  compatible  with  the 
data  of  McKague  [40]  and  satisfy  the  theoretical  relation 
derived  in  Ref.  [41,  p.  69]. 

8 . 3  Modeling  of  Epoxy  Properties 


The  observations  of  the  preceding  section  are  used 
for  modeling  the  material  properties  of  epoxy  that  are 
given  by  Eq .  (8.1). 

The  glass  transition  region  of  epoxy  resin  is  not 
broad  [6],  therefore,  a  glass  transition  temperature  is 
used  instead.  The  temperature  T  is  usually  obtained  by 
measuring  the  expansion  of  a  specimen  as  function  of 
temperature.  The  point  where  the  epoxy  stops  expanding  as 
temperature   increases  corresponds   to   the  first  deviation 

from  the  glassy  state  and  is  termed  T  . 

g 

According  to  the  experimental  data  plotted  in 
Figure  8.4,  T  is  stongly  dependent  on  absorbed  moisture. 
These  results  show  that,  as  the  moisture  content  of  epoxy 
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increases,   the  transition  temperature  moves  to  the  left  in 

Figures  8.1-8.3.  Hence,   the  abrupt  change  of  the  material 

properties  starts  at  a  lower   temperature  as   the  moisture 

content  increases.  This  fact  and  the  conclusions  reached  by 

previous   investigators  [12]   suggest   that   the   following 

modelings  of   E' ,  u ' ,  and  n '   are  appropriate, 
to        mm        m 


E'  =  E"  f 
m     o 


T  -  T 


T  -  T 
g    o 


(8.2) 


v       =    v       g 
m     o 


T  -  T  ■» 


T  -  T 
g    o 


(8.3) 


■n '  =  n '  h 
m     o 


T  -  T 


T  -  T 
g    o 


(8.4) 


where   the   temperature  T    is  equal   to  273K.   The  moisture 

o 

concentration   appears    implicitly    in    T  .    The    glass 
transition  temperature  is  represented  by 


T   =  210  exp(-  9c)   (  C) 


(8.5) 


where  c  is  the  moisture  concentration. 

This  modeling  has  been  chosen  so  that  it  does  not 
represent  the  material  properties  beyond  T  ,  since  the 
study  of  epoxy  in  the  rubbery  stage  is  not  within  the  scope 
of  this  research.  Equations  (8.2)-(8.4)  are  valid  only  for 
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the    continuous    parts    of    the    curves    plotted    in 
Figures  8.1-8.3.  , 

8 . 4  Resu 1 t s 

All  test  specimens  are  conditioned  in  a  constant 
relative  humidity  environment  until  moisture  equilibrium  is 
reached.  Then,  the  test  pieces  undergo  the  impulse  hammer 
technique  and  the  four-point  flexure  tests  to  determine  the 
storage  moduli,  the  material  damping,  and  Poisson's  ratio 
at  several  temperature  and  moisture  contents. 

8.4.1  Complex  Moduli  of  Epoxy 

Storage  modulus.  The  experimental  data  on  the  storage 
modulus  of  epoxy  in  term  of  temperature  at  three  different 
equilibrium  moisture  concentrations  are  plotted  in 
Figure  8.5.  It  can  be  concluded  that  increase  in  either 
temperature  or  moisture  content  or  both  results  in  a 
decrease  in  the  storage  modulus.  Plotting  these  data  in 
terms  of  moisture  content  in  Figure  8.6  does  not  lead  to 
any  additional  insight.  But,  representing  these  results  in 
term  of  the  following  normalized  non-dimensional 
temperature 


T  -  T 

T   =  ! 

n    T  -  T 

g    < 


(8.5) 
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in   Figure  8.7   shows   a   clear   trend.   Experimental   studies 

have  shown  that  the  modulus  of  polymer  is  very  low  at  the 

glass  transition  temperature,  therefore,  adding  the  value 

E   =0  for  T  =  T   to  the  data  yields  the  following  modeling 
m  g  &         & 


E1  =  4.0(1  -  T  ) 
m       v      n' 


(GPa) 


(8.6) 


Material  damping.  Similarly,  the  experimental  data  of 
the  hygrothermal  effects  on  the  damping  of  epoxy  are 
plotted  in  three  Figures  (8.8-8.10).  There  is  very  little 
change  in  damping  for  all  the  considered  conditions. 
Therefore,  it  is  proposed  to  let 


tj   =  0.018 
m 


(8.7) 


for  temperatures  up  to  80  C  and  moisture  contents  up  to  5%. 
The  conclusion  that  the  hygrothermal  effects  on  the  damping 
of  epoxy  is  negligible  is  qualitatively  corroborated  by 
Putter  et  al  .  [38].  A  quantitative  comparison  cannot  be 
made  since  these  researchers  have  not  included  in  their 
publication  the  values  of  the  fiber  volume  fraction  and 
moisture  content  of  the  test  specimens. 


Po  i  s  son ' s   ra  t  i  o .   The   experimental   values   of   the 
Poisson's  ratio   in  terms  of   temperature  at   two  different 
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moisture  contents.  are  plotted  in  Figure  8.11.  These 
results  show  that  Poisson's  ratio  increases  at  a  negligible 
rate  as  temperature  varies  from  0  to  80  C.  Representing  the 
same  data  in  terms  of  the  moisture  content  up  to  M  =  4.5% 
(Figure  8.12)  shows  that  the  moisture  effect  is  also 
negligible.  Therefore, 


v'     =0.32 
m 


(8.8) 


for  temperatures  up  to  80  C  and  moisture  contents  up  to  5%. 

Since    v'         equals  0.5  at  the  glass  transition  temperature 
m 

(T   =  0),  the  plot  of  Poisson's  ratio  versus  the  normalized 
v  n     ' 

temperature  has  been  extrapolated  as  shown  in  Figure  8.13. 
The  extrapolation  displays  a  qualitative  trend  only. 

8.4.2  Complex  Moduli  of  Composites 


The  complex  moduli  of  Glass/Epoxy  and  Graphi te/Epoxy 
in  terms  of  moisture  content  and  temperature  can  be 
determined  by  using  the  fibers'  properties  given  in 
Table  7.2,  Eqs .  (8.6)  through  (8.8)  and  the  mi cr omechanics 
formulas  (Eqs.  (3.20)). 

This  procedure  is  illustrated  by  determining  the 
storage  moduli  and  the  damping  of  a  Glass/Epoxy  lamina  with 
a  fiber  volume  fraction  of  0.5  and  a  Graphi te/Epoxy  lamina 
with  a  fiber  volume  fraction  of  0.7. 
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G  las  s/Epoxy  .  The  parameters  E  '  .  ,  E'_,  G '  .  rj1  -  , 
^9  9'  ^r  '  anc*  u  1  9  versus  the  normalized  temperature  are 
plotted  in  Figures  8.14-17.  The  experimental  data 
substantiate  the  theoretical  results. 

Graph  i  te/Epoxy.  Similarly,  E  '  ..  ,  E'„,  G  -'  _ ,  1 1  i  ■ 
^99'  ^r  •  an<i  u  i  9  versus  the  normalized  temperature  of 
Graphi te/Epoxy  are  plotted  in  Figures  8.18-21. 

For  both  Glass/Epoxy  and  Graphi te/Epoxy ,  the  results 
show  that  the  matrix-dominated  parameters  (E'  and  G'  )  are 
strongly  affected  by  moisture  and  temperature,  while  the 
fiber-dominated  parameters  (E'  ,  u'  )  stay  practically 
cons  tan  t . 
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Fig.  8.1    Schematic  variation  of  the  storage  modulus  of 
epoxy  with  temperature. 
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Fig.  8.2    Schematic  variation  of  Poisson's  ratio  of  epoxy 
with  temperature. 
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Fig.  8.3    Schematic  variation  of   damping  of   epoxy  with 
temperature . 
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Fig.  8.4    Glass   transition   temperature   of   epoxy 
Delasi  and  Whiteside  [6]. 
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Fig.  8.5  Experimental  data  of  the  storage  modulus  of 
epoxy  as  a  function  of  temperature  at  diverse 
constant  moisture  contents. 
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Fig.  8.6  Experimental  data  of  the  storage  modulus  of 
epoxy  as  a  function  of  moisture  content  at 
diverse  constant  temperatures. 


70 


*     Experimental   data 


Fit  to  data 


0.2        0.4         0.6         0.8  1 

Normalized  temDerature 


Fig.  8.7    Experimental   data   of   the   storage   modulus   of 
epoxy  as  a  function  of  normalized  temperature 
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Fig.  8.8  Experimental  data  of  damping  of  epoxy  as  a 
function  of  temperature  at  diverse  constant 
moisture  contents. 


0.03 


oi   0.025  \- 
_c 
5. 

E 
o 

a 


0.02 


0.01 5  C  J- 


0.01     - 


0.005 


a 


'■■■'■ 


*  T  =  20  C 
O  T  =  50  °C 
Q     T  -  70  °C 


2  3  4  5 

Moisture  content  {%) 


Fig.  8.9  Experimental  data  of  damping  of  epoxy  as  a 
function  of  moisture  content  at  diverse 
constant  temperatures. 
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Fig.  8.10   Experimental   data  of   the   storage   modulus   of 

epoxy  as  a  function  of  normalized  temperature 
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Fig.  8.12   Experimental  data  of  Poisson's  ratio  of  epoxy 
in  term  of  moisture  content. 
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Fig.  8.17   Poisson's  ratio    ( v ' „ )    of  Glass/Epoxy  versus 


T   =  (T  -  T  )/(T   -  T  ) . 
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CHAPTER  9 
HYGROTHERMAL  EFFECTS  ON  STRESS  FIELD 


9 .  1  In  t  roduc  t  i  on 

The  hygro thermal  effects  on  the  stress  field  are 
investigated  by  considering  an  infinitely  long,  finite 
width  and  symmetric  composite  laminate  undergoing 
hygrothermal  loadings.  The  Finite  Element  Method  is  used  in 
order  to  estimate  the  magnitude  of  hygrothermal  stresses  in 
laminated  composites  (see  Appendix  B).  The  geometry  of  a 
laminate  and  the  finite  mesh  of  a  quarter  cross-section  are 
shown  in  Figure  9.1  and  the  boundary  conditions  are  given 
by 

v  =  0    for  (y , z)  =  (0,  z) 

(9.1) 

w  =  0    for  (y , z)  =  (y  ,0) 

where  v  and  w  are  the  displacements  in  the  y  and  z 
directions,  respectively.  The  grid  consists  of  24  eight 
node  isoparametric  elements  and  93  nodes.  Only  24  elements 


79 


80 

are  used  since  increasing  the  number  of  elements  to  48 
results  in  a  relatively  small  change  in  the  stress 
magnitudes.  The  material  properties  in  terms  of  temperature 
and  moisture  content  have  been  derived  in  the  preceding 
chapter.  The  constitutive  equations  are  given  by  Eq .  (B.12) 
and  can  be  written  in  matrix  form  as 

{a}    =  [Q]({€>  -  {a}AT  -  (P}c)  (9.2) 

where  {a}  and  {(3}  are  the  vectors  of  thermal  and  moisture 
expansion  coefficients. 

9 . 2  Description  of  Study  Cases 

The  considered  stacking  sequence  is  the  [(90/0)9]„ 
lay-up.  The  cross-ply  laminate  is  preferred  over  other 
laminate  since  hygrothermal  loadings  induce  very  high 
stresses  in  this  case.  The  volume  fiber  fractions  of  the 
Glass/epoxy  and  the  Graphi te/Epoxy  are  0.5  and  0.7, 
respectively.  The  thickness  and  the  width  of  the  laminates 
are  assumed  to  be  2  mm  and  20  mm,  respectively. 

Three  cases  of  moisture  gradients  are  applied.  They 
are  described  in  Figure  9.2  and  Table  9.1.  Cases  A  and  C 
correspond  to  the  dry  and  moisture  saturated  states, 
respectively.  While  the  non-uniform  moisture  gradient 
(case  B)  corresponds  to  a  moisture  profile  as  derived  in 
section  2.3.  Two  uniform  temperatures  (20  C  and  80  C)  are 
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used.  All  laminates  are  assumed  to  be  initially  (dry  at 
20  C)  free  of  stress.  Hence,  residual  stresses  are  not 
taken  into  account.  The  elastic  moduli  used  in  computing 
the  stresses  are  approximated  by  the  real  parts  of  the 
complex  moduli.  Therefore,  the  hygrothermal  effects  on  the 
elastic  properties  can  be  deduced  from  the  results  given 
in  Chapter  8 

9 . 3  Numerical  Results  and  Discussion 

For  all  considered  cases,  the  following  remarks  can 
be  drawn:  at  z/h  =  constant,  the  stresses  away  from  the 
free  edge  stay  constant  and  the  shear  stress  (a  )  is 
zero,   but,   as  y/b  approaches  1,    a  takes  significant 

non-zero   values   and   there   are   small   variations   in   the 

values  of  the  other  stresses.  Hence,   the  stresses   a     ,  a 

y    z 

and  a  are  plotted  across  the  section  of  the  laminate  at 

x 

y/b  =  0.472  and  the  shear  stress   a  is  plotted  across 

yz 

the  section  at  y/b  =  0.993  (close  to  the  free  edge). 

The   stresses  are  compared  to   typical  strengths  of 

Glass/Epoxy  and  Graphi te/Epoxy  that  are  provided  in 
Table  9.2. 

9.3.1  G lass/Epoxy 


The  equilibrium  moisture  concentration,   c  ,   of  the 
Glass/Epoxy  material  is  0.025. 
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The  stress  a  is  plotted  in  Figure  9.3.  It  reaches 

y 

a  maximum  magnitude  of  166  MPa.   for  case  C  at  20  C.   It  is 

compressive  for  the  0   layer  and  tensile  for  the  90   layer. 

The   stress     a  is   shown   in   Figure  9.4.   It   is 

compressive  everywhere  and  reach  a  magnitude  of  288  MPa. 

for  the  case  C  at  20  C.  The  stress  a       is  also  compressive 

x 

(Figure  9.5)  and  reaches  a  maximum  of  245  MPa.  .  The  free 
edge  shear  stress    a  (Figure  9.6)  is  very  significant 

since  its  maximum  magnitude  is  about  80  MPa. . 


9.3.2  Graphi  te/Epoxy 

The  equilibrium  moisture  concentration   c    for  these 

^  CO 

cases  is  0.015.   The  stresses    o    ,  a    ,  a    ,   and  a  are 

y    z    x         yz 

plotted  in  Figures  9.7-10.  These  results  show  the  same 
trend  as  for  the  Glass/Epoxy  cases.  However,  since  the 
moisture  concentration  is  lower  and  graphite  fibers  have 
stiffer  moduli  and  lower  coefficient  of  thermal  expansion, 
the  magnitude  of  the  stresses  is  smaller. 

9.3.3  Summary 

The  hygrothermal  conditions  used  in  the  preceding 
sections  are  practically  achieved  only  under  very  adverse 
conditions.  Hence,   the  induced  stresses  can  be  considered 
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an  upper  bound  for  hygro therma 1  stresses.  The  results  yield 
the  following  observations: 

1)  The  stresses  induced  by  temperature  only  (dry  at 
80  C)  are  much  smaller  than  those  induced  by  high  moisture 
content . 

2)  The  stresses  due  to  a  non-uniform  moisture 
gradient  can  be  as  high  as  those  induced  by  the  saturated 
mo  isture  case . 

3)  Since  the  hygrothermal  conditions  degrade  the 
modulus  of  the  epoxy  matrix,  the  stresses  caused  by  the 
most  severe  hygrothermal  condition  (moisture  case  C  at 
80  C)  are  lower  than  for  some  of  the  other  cases. 

4)  The  hygrothermal  stresses  of  the  cross-ply 
laminates  are  very  significant  since  their  magnitude  is  of 
the  same  order  of  those  of  the  strengths  given  in 
Table  9.2. 


S4 


Table  9.1   Description  of  cases  in  Figure  9.2 


Case 


Description 


A  Dry 

B      (M  -  M.)/(Mro-  M.)  =  0.5  (absorption  cycle) 

C      (M  -  M.)/(Mm-  M.)  =  1.0  (fully  saturated) 


Table  9.2   Typical     strengths 
Graph i  te/Epoxy 


of     Glass/Epoxy 


and 


Strengths 
(MPa. ) 


G 1  as  s/Epoxy 


Graphi  te/Epoxy 


1000.0 
1000.0 

30.0 
140.0 

40.0 


1200.0 

700.0 

40.0 

70.0 

70.0 


X   (X  )  =  Longitudinal  strength  in  tension  (compression) 
Y   (Y  )  =  Transverse  strength  in  tension  (compression) 
S  =   Shear  strength. 
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Fig.  9.1    Geometry  of  a  laminate  and  finite  mesh  of  a  1/4 
cross-section  area. 
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Moisture  case    A 


Moisture  case  B 
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Fig.  9.2    Description  of  the  applied  moisture  gradients 
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Fig.  9.10   Profile  of  the  hygro thermal  stress  a  across 

a    [(90/0)_]     Graphi te/Epoxy    laminate    at 
y/b  =  0.993. 


CHAPTER  10 
HYGROTHERMAL  EFFECTS  ON  COMPLEX  STIFFNESSES 


10.1  Introduction 


The  hygrothermal  effects  on  the  in-plane  (A.  .)  and 

the  bending  (D.  .)  complex  stiffnesses  of  Glass/Epoxy  and 

Graphi te/Epoxy    angle-ply      [(iG)„]„       laminates    are 

investigated.  The  applied  moisture  gradients  are  the  cases 

A,  B,  and  C  that  are  given  in  section  9.2  and  the  uniform 

applied     temperatures   are   20  C   and   80  C.   The   material 

properties   in   terms   of   moisture  content  and   temperature 

have  been  determined  in  Chapter  8.   The  fibers  properties 

are  given  in  Table  7.2.  The  theoretical  expressions  of  A.  . 

1J 

and  D.  .  in  terms  of  the  fibers  and  matrix  properties  have 
been  developped  in  Chapter  5  and  Appendix  A. 


10.2  Numerical  Results  and  Discussion 


10.2.1  Glass/Epoxy 


The  fiber  volume  fraction  of  the  Glass/Epoxy  is  0.5 
and  the  equilibrium  moisture  content    cro    is  0.025.  Since 
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the   fibers   do   not   absorb   any   moisture,   the   equilibrium 

moisture  concentration  of  the  matrix  is  0.05.  The  thickness 

of   the  laminate  is  2.0  mm  and  A.  .  and  D.  .  are  normalized 

i  J         ij 

with    respect    to    75.86  x  10     N/m    and    25.29    N.m, 
respec  t  i ve ly . 

The  real  part  of  the  longitudinal  in-plane  complex 
stiffness,  A',  and  its  corresponding  damping ,  .i)  are 
plotted  in  Figure  10.2.  As  the  moisture  gradients  and  the 
temperature  change,  the  relative  changes  of  A'  vary  from 
6%  (for  ±9  =  0°)  to  33%  (for  ±9  =  90°).  The  line  style 
legend  of  the  figures  of  this  chapter  is  defined  in 
Figure  10.1. 


Similarly,  A'   and   .?).„   are  plotted  in  Figure  10.3 

and  A'   and    tticc         in  Figure  10.4.  All  these  cases  show 
oo         1  bo 

that  the  hygrothermal  effects  on  A.  .  is  matrix  dominated. 

The  real  part  of  the  bending  stiffnesses,  D.  .,  and 
their  corresponding  damping,  -n .  .,  are  plotted  in  Figures 
10.5-7.   These  results  yield  similar  conclusions   to   those 


of  A 


ij 


10.2.2  Graphi  te/Epoxy 


The  volume  fraction  of  the  Graphi te/Epoxy  laminate  is 
0.7  and   its   equilibrium  moisture  concentration    c     is 

^  CO 

0.015.  Since  all  the  moisture  is  absorbed  by  the  matrix  , 
cro   of  the  epoxy  is  0.050.  The  laminate  is  2.0  mm  thick. 
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The  real  parts  of   the  complex  stiffnesses  (A!  .  and 

ij 

D!  .)  and  their  corresponding  damping   (TT7.  .  and  „n.  .)   are 

plotted   in   Figures   10.8-13.   The   terms   A!  .   and   D!  .   are 

ij         ij 
/? 

normalized  with  respect  to  311.3  x  10  N/m  and  103.77  N.m, 
respectively.  These  results  show  the  same  tendency  as  those 
of  Glass/Epoxy.  Since  the  volume  fraction  and  the 
longitudinal  modulus  of  the  graphite  fibers  are  higher,  the 
hygrothermal  effects  are  less  pronounced. 


10.2.3  Summary 


The  effects  of  moisture  on  the  stiffnesses  A. .'s  and 

i  J 

D. .'s  of  composites  at  room  temperature  are  negligible  for 
all  the  considered  cases.  But,  as  temperature  increases, 
the  combined  influence  induces  significant  changes  in  the 
complex  stiffnesses  especially  for  the  matrix 
dominated  terms. 
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Moisture  gradient  case  A  at  20°C 

Moisture  gradient  case  B  at  20°C 

Moisture  gradient  case  C  at  20°C 

Moisture  gradient  case  A  at  80°C 


Moisture  gradient  case     B    at  80°C 


Moisture  gradient  case     C    at  80  C 


Fig.  10.1    Line  style  legend  of  Figures  10.2-13 
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Fig.  10.2   Complex  in-plane  stiffness  A11  of  Glass/Epoxy. 

a)  Non-dimensional  Real  part  ;  b)  corresponding 
damping . 
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Fig.  10.3   Complex  in-plane  stiffness  A19  of  Glass/Epoxy. 

a)  Non-dimensional  Real  part;   b)  corresponding 
damping . 
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a)  Non-dimensional  Real  part;   b)  corresponding 
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Fig.  10.6   Complex  bending  stiffness  D19  of  Glass/Epoxy. 

a)  Non-dimensional  Real  part;   b)  corresponding 
damping . 
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a)  Non-dimensional  Real  part;   b)  corresponding 
damping . 
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CHAPTER  11 
CONCLUSION 


Theoretical  and  experimental  methods  have  been 
incorporated  in  order  to  determine  the  effects  of 
temperature  and  absorbed  moisture  on  the  complex  moduli  of 
composite  materials.  The  effects  of  hygrothermal  loadings 
and  of  the  changes  of  the  complex  moduli  on  the  stress 
field  and  on  the  structural  damping  of  composite  laminates 
are  also  analyzed. 

The  Fickian  theory  of  mass  diffusion  has  been  used 
for  analyzing  the  diffusion  of  moisture.  The  Fickian 
diffusion  is  adequate  for  the  experimental  determination  of 
moisture  diffusion  through  unstressed  test  specimens. 
Therefore,  theories  that  incorporate  the  coupling  of 
moisture  diffusion  with  stress,  viscoelastic  relaxation, 
entropy  inequality,  etc.  [43-45],  have  not  been  included. 

The  complex  moduli  of  unidirectional  composites  are 
expressed  in  terms  of  the  constituent  material  properties 
in  Chapter  3  by  using  the  following  steps: 

-  first,  the  elastic  moduli  of  unidirectional 
composites  in  terms  of  the  fibers  and  epoxy  properties  are 
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obtained  by  using  the  rule  of  mixture  and  the  Halpin-Tsai 
equa  t  i  ons , 

-  then,  the  correspondence  principle  is  applied  to 
determine  the  complex  moduli. 

Equations  (3.20)  show  that  in  order  to  determine  the 
hygrothermal  effects  on  the  complex  moduli  of  composites, 
the  effects  on  eight  distinct  parameters  need  to  be 
assessed.  But,   these  effects  on  the  fibers  properties  are 

negligible.  Hence,  only  three  terms,    E'(T,c),  t?  (T,c)  and 

&   &  J  mv     '  mv    ' 

u'(T,c),   have  to  be  experimentally  measured.  The  necessary 

test  procedures  are  descibed  in  Chapter  6. 

The  complex  moduli  of  epoxy  in  terms  of   temperature 

and   moisture   content   are   presented   in   Chapter  8.   The 

storage  moduli,   E',   is   strongly  dependent  on   temperature 

and  moisture  content.  But,   n   and  v'        stay  constant  up  to 

m        m 

a  moisture  content  M  =  4.5%  and  a  temperature  T  =  80  C. 
More  severe  hygrothermal  conditions  could  not  be  reached  in 
the  environmental  chamber  due  to  operating  temperature 
limitations  of  the  chamber  and  of  the  motion  and  force 
transducers  used  in  the  impulse  hammer  technique. 

The  experimental  results  of  the  epoxy  properties  in 
combination  with  the  mi c romechani cs  formulas  of  Chapter  3 
are  used  to  determine  the  hygrothermal  effects  on  the 
complex  moduli  of  unidirectional  Glass/Epoxy  and 
Graphi te/Epoxy  laminates.  It  is  shown  that  only  the  matrix 
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dominated  terms  (E99  and  G'  )  are  strongly  affected  by 
temperature  and/or  moisture. 

Hygrothermal  conditions  and  changes  in  the  complex 
moduli  influence  the  stress  field  and  the  structural 
damping  of  laminated  composites. 

The  stresses  induced  by  hygrothermal  conditions  in 
Glass/Epoxy  and  Graphi te/Epoxy  cross-ply  laminates  are 
illustrated  in  Chapter  9.  The  [(90/0)9]~  lay-ups  are  chosen 
since  they  display  the  highest  hygroscopic  stresses.  The 
induced  stresses  in  the  90  layers  are  of  the  same  order  of 
magnitude  as  the  strengths.  Hence,  high  loadings  induced  by 
temperature  and  moisture  content  can  lead  to  failure. 

The  complex  stiffnesses  and  structural  damping  of 
[(±6)„]~  Glass/Epoxy  and  Graphi te/Epoxy  laminates  are 
analyzed  in  Chapter  10.  These  results  display  the  same 
general  trend  as  those  of  the  complex  moduli  of 
unidirectional  composites.  That  is,  only  the  matrix 
dominated  terms  are  strongly  affected  by  moisture  and/or 
temperature . 

The  important  conclusions  of  this  investigations  are: 

•  hygrothermal  effects  are  very  significant  for  the 
matrix  dominated  properties  only; 

•  for  certain  laminates,  severe  hygrothermal  stresses 
alone  can  lead  to  failure.  Therefore,  hygrothermal  effects 
need  to  be  taken  into  account  when  designing  composite 
structures  that  are  subjected  to  moisture  and/or 
tempera  ture . 
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Recommendations  and  additional  remarks  that  arose 
from  this  investigation  are  listed  below. 

•  The  results  and  conclusions  could  vary  among 
different  material  systems.  Hence,  it  might  be  necessary  to 
repeat  the  tests  and  the  methodology  for  different 
mater  ial s . 

•  The  experimental  results  are  valid  up  to  a  80  C 
temperature  and  a  4.5%  moisture  content.  In  order  to 
determine  the  properties  of  epoxy  through  a  wider  range  of 
temperature  (from  -  50  C  to  200  C)  and  moisture  content,  an 
environmental  chamber  and  transducers  that  can  operate 
under  more  severe  hygrothermal  conditions  need  to  be  used. 


APPENDIX  A 
COMPLEX  STIFFNESSES  OF  COMPOSITES 


A . 1  Elastic  Stiffnesses 

The  in-plane,  coupling  and  bending  stiffnesses  of  a 
general  laminated  elastic  composite  plate  (Figure  A.l)  are 
given  by 


A.  .  = 


•+h/2 


Q.  .  dz 


-h/2 


r+h/2 


B 


i  J 


zQ.  .  dz 
i  J 


-h/2 


(A.l) 


D. 


ij 


■+h/2 


-h/2 


z2Q.  .  dz 


where  the  components  of  the  transformed  matrix 


[Q]  = 


^11   ^12   ^16 
Q12  Q22  Q26 

Q16   Q26   Q66 


(A. 2) 
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are    given    by 


Q         =    Qnm4    +    2(Q12    f    2Q66)m2n2    ♦    Q^n4 


Q12    =    <Q11    +    Q22    "    4Q66)m2n2    +    Q12^4    +    n^ 


2    2 


Q22    =    Qlln       +    2^Q12    +    2(Wm    n       +    Q 


22 


m 


(A. 3) 


Q16    =     (Qn    -    Q12)m3n    +     (Q12    -    Q22)mn3    +    2Q66(m2-    n2)mn 


3  3  2  2 

Q26    =    ^Qll    "    Ql2^mn      +    ^Q12    "    Q22^m    n    +    2Q66^m       ~    n    ^mn 


Q66    =    <Q11    +    Q22    "    2Q12    "    2Q66)m2n2    +    Q66^4    "    ^ 


and 


m    =    cos    0 


n    =    sin    0 


The  angle    0    represents   the  fibers  orientation  of   the 
lamina  under  consideration. 


The  parameters  Q.  .'s  can  be  expressed  in  terms  of  the 
longitudinal  (E11)  and  transverse  (EQO)  Young  modulus,  the 


11 


22 


shear  modulus  (G1<?)  and  the  major  Poisson's   ratio  ( v  1 ~ ) 
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Qll  =  Ell/(1  "  u12u21> 


Q12  =  D12E11/(1  "  U12U2l) 


Q22  =  E22/(1  "  U12U2l) 


(A. 4) 


Q66  "  G12 


U21  =  U12E22/E11 


The  parameters  E1  1  ,  E99  •  G.„  and  uio  are  obtained  as 
functions  of  the  properties  of  the  constituent  materials 
(fibers  and  matrix)  by  using  Eqs .  ( 3 . 5) - (3 . 10) . 


A . 2  Complex  Stiffnesses 


The  laminate  complex  stiffnesses  are  determined  by 
carrying  the  following  steps 

i)     The     e las t i c-v i scoe las t i c     correspondence 
principle  is  applied  to  Eqs.  (3.5)-(3.10) 

ii)   The   resulting   complex   values   replace   their 
corresponding  elastic  modulus  in  Eqs.  (A. 4) 

iii)  Consequently,  the  complex  stiffnesses  Q.  .'s  are 
obtained  and  substituted  for  their  correspondig  elastic 
stiffnesses  in  Eqs.  (A. 3) 
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iv)   Finally,   the   complex   transformed   stiffnesses 

Q.  .'s  are  substituted  in  Eqs .  (5.1). 
i  J 

This  procedure  is  easily  executed  with  a   FORTRAN  program. 


APPENDIX  B 
DEVELOPMENT  OF  THE  FINITE  ELEMENT  METHOD 


The  finite  element  displacement  method  used  to  derive 
the  stress  field  in  the  cases  of  Chapter  9  is  presented  in 
this  Appendix.  The  F.E.M  presented  below  is  used  only  to 
give  a  first  approximation  of  the  hygrothermal  stress 
magnitudes  in  laminated  composites.  In  order  to  obtain  more 
accurate  results,  another  method  that  takes  into  account 
the  zero-stress  boundary  conditions  on  the  free  surfaces 
should  be  applied. 

B  .  1  Equilibrium  Equations 

The  stress-strain  relation  of  a  linear  elastic  solid 
continuum  undergoing  any  type  of  loading  can  be  written  in 
the  form 


{a}    =  [D]({€}  -  {€°})  +  {a°}  (B.l) 


where  [D]  is  the  elasticity  matrix,  {€}  is  the  strain 
vector,  {€  }  is  the  hygrothermal  and/or  initial  strain 
vector,   {a    }   is  any  initial  stress  vector. 
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In  this  method,  it  is  assumed  that  the  displacements 
have  unknown  values  only  at  the  nodal  points.  The 
displacements  are 


n 


{6}    =    [N]{66}  =  I    [N.]{6.} 


(B.2) 


i  =  l 


where 


[N]  =  [Nr  N2.    ,Nn] 


[N.]  =  N.[I] 


and  n  is  the  number  of  variables  per  node  and  [I]  is  the 
identity  matrix. 

The  strains  in  terms  of  the  displacements  are  given 
by  the  following  expression 


{€}  = 


■ 

€ 

y 

e 

z 

av 

•     =  ■ 

ay 

3w 
dz 

7 

yz 

dv 
.    3z 

n 


=  [B]{66}  =  I    [B.]{6.} 


i  =  l 


(B.3) 


where  the  strain  matrix  is  defined  as 


[B.]  = 
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3N. 

i_ 

ay 
o 

9N. 


0 

3N. 

i_ 

dz 

3N. 

i 


3z     5y 


(B.4) 


The  set  of  functions  N.  are  called  the  shape  functions  and 
are  subsequently  defined. 

The  continuum  is  subdivided  into  a  finite  number  of 
elements.  Consider  an  element  that  is  acted  upon  by  nodal 
forces  {F  }  and  body  forces  {p}.  Then,  application  of  the 
virtual  work  principle  to  an  element  e  yields 


[5e]T{Fe}  + 


[6]'{p}dV  = 


[€]'{CT}dV       (B.5) 


Substitution  of  Eqs.  (B.1)-(B.4)  into  Eq .  (B.5)  results  in 


{Fe}  + 


[N]1{p}dV  - 


[B]*[D][B]dV 


{5e} 


[B]T[D]{€°}dV  + 


[B]T{a°}dV    (B.6) 


or 


{Fe}    +    {Fp    +    {F^o}    +    {FeaQ}    =    [Ke]{5e}  (B.7a) 


where 
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[Ke]  = 


T 
[B]  [D][B]dV   =  element  stffness  matrix 


{Fe}  = 
P 


[N]1{p}dV 


=  equivalent  nodal  body  force 


V 


IF   >  = 


[B]  [°]{e  )dV  =  hygrothermal  or  initial 
e  strain  loading 


<FCTO}  =  " 


T 
[B]  {a}dV    =  initial  stress  loading 


These  equations  are  valid  for  one  element  only.  For  the 
complete  structure,  Eq .  (B.6)  should  be  summed  over  all 
elements.  It  is  noted  that  if  an  element  is  subjected  to 
surface  traction  forces,  {t},  then  the  additional 
equivalent  nodal  force 


{Fp  = 


[Nr{t}dS 


(B.7b) 


should  be  added  to  the  left  hand  side  of  Eq .  (B.7a) 


B . 2  Program  Organization 


Hinton  and  Owen  [42]  have  developed  a  detailed  F.E.M. 
code  to  solve  isotropic  beam,  plane  stress/strain  and  plate 
bending  problems.  Their  procedure  for  the  plane 
stress/strain  has  been  modified  so  that   it   is  adapted  to 
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composite  laminates  that  have  nonconstant  material 
properties  and  undergo  thermal  and  moisture  induced 
expans  i  on . 

The  F.E.M.  operations  are  performed  by  modular 
subroutines.  The  general  organization  of  these  programs  is 
shown  in  Figure  B.l. 

B . 3  Shape  Functions,  Jacobian  and  Strain  Matrix 

Eight  node  isoparametric  elements  are  used  in  the 
F.E.M.  code.  In  this  case,  the  shape  functions  are  used  to 
approximate  both  the  geometry  and  the  displacement  field. 
The  coordinates  of  any  point  of  the  element  shown  in 
Figure  B . 2  are 


y(f.T])  =  J   N.(f,n)y. 


i  =  l 


8 


(B.8) 


ZU.TJ)   =  I     N.(f,Tj)z. 

i  =  l 


where  (y.,  z.)  are  the  coordinates  of  the  node  i  and  the 
quadratic  shape  functions  are  defined  as 


Nj  =  -  (1  -  f)(l  -  77)  (  1  +  I    +    T7)/4 
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n0  =  (i  -f  )(i  -  n)/2 


N3  =  (1  +  f)(l  -  i7)(f  -  t!  -  l)/4 


N4  =  (1  +  f)(l  -  T|  }/2 


(B.9) 


If-  =  (1  +  f)(l  +  Tj)(f  +  n  -  l)/4 


N6  =  (1  -  f  )(1  +  i?)/2 


N?  =  (1  -  f)(l  +  tj)(-  f  +  T7  -  l)/4 


N8  =  (l  -  HO  -  n  )/2 


The  direction  of  the  local  curvilinear  coordinates  f  and  n 
are  given  in  Figure  B.2.  The  displacements  at  any  point  are 
expressed  as 


8 


v(f.T7)  =  I    N.(f.T,)vi 


i  =  l 


8 


(B.10) 


w(f.n)  =  J    N.(f,77)w. 


i  =  l 


The  Jacobian  matrix  [J(f .l)]  is  necessary  to  derive 
the  elemental  volume  dV ,  and  area  dS .  It  is  defined  as 
f o 1  lows 
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[J]  = 


9y 

ay 
377 


dz 

as 
dn 


i  =  l 


3N. 

aF^i 


aN. 
i 

dr)  yi 


9N. 

arzi 


aN. 
i 

a z  • 

OT)         1 


(B.ll) 


Then  the  elemental  area  dS  is  given  by 


dS  =  dy  dz  =  det[J]  df  dn 


Once   the   shape   functions  are  chosen.   the   strains  can  be 
written  in  terms  of  N.  and   6.   (Eq.  (B.3)). 


B . 4  Elasticity  Matrix 

The  stress-strain  relation  of  an  orthotropic  lamina 
is  g  i ven  by 
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(B.12) 

where   the   transformed   stiffnesses   that   have   not   been 

defined  in  Appendix  A  are  expressed  as 
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Q13  =  Q13m2  +  Q23r2 


Q23  =  Q13n2  +  Q23m2 


Q33  "  Q33 


(B.13) 


Q36  =  (Q13  "  Q23)mn 


Q44  =  G^m2  +  G12n2 


Q45  =  (G12  "  G23)mn 


Q55  =  G12m   +  G23n 


The  stiffnesses  Q.  .'s  are 


Q13  =  U21^  +  U12>E11/Q 


Q23  =  U12^  +  U21^E22/Q 


Q 


33 


(1  -  vl2v2l)E22/Q 


Q  = 


=  l    -    2l,12u21^  +  2u23> 


-  i) 


23 


where 
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r    v 
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3  -  4u   +  G  /Gnr, 
m m   12 
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U23  =  Um 


For  the  case  of  Chapter  9,  Eq .  (B.12)  can  be  reduced  to 
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yzj 
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yzj 
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x       x 
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y     y 

a  AT  +  6  c 
z       z 

a   AT  +  0   c 
L  xy      ^xy 


(B.15) 


{a}    =    [D1]{€}  -  [D2]{C°} 


The  last  equation  is  a  reduced  form  of  Eq .  (B.l)  with  the 
initial  stresses  left  out. 

The  through  the  thickness  stress  is  given  by 


a      =    Qni(-  a    AT  -  B    c)  +  Q10(€   -  a    AT  -  B    c) 
x     llv    x      'x-'     12vy     y       y 


+  Q.,(€^  -  a  AT  -  Be)  +  Q1R(-  ov„AT  -  Be)   (B.16) 


13   z     z 


16v    xy 


xy 


12S 


B . 5  Element  Stiffness  Matrix 

A   submatrix   of   the   element   stiffness   matrix   [K 
linking  nodes  i  and  j  is  evaluated  from  the  expression 


[K  ,  .]  = 


■+1 


■+1 


-1J-1 


[B.]  [D][B.]  t  det[J]  df  dT]    (B.  17) 


where  t  is  the  thickness  of  the  element  under 
consideration.  The  integration  is  done  by  using  a  3-point 
Gauss  integration  rule. 

B . 6  Equivalent  Nodal  Loadings 

B . 6 . 1  Element  Edge  loadings 

An  element  edge  might  have  both  tangential  and  normal 
distributed  load  per  unit  length.  Every  edge  of  the 
isoparametric  element  has  three  nodal  points.  The  values  of 
the  normal  and  tangential  loads  at  each  nodal  point  are 
called  (p  ).  and  (p  )..  Then,  the  distributed  loads  at  any 
point  along  the  edge  are  given  by 


pn 


"  I'i 


i  =  l 


(Pn>i 


<Pn>i 


(B.18) 
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It   can   be   shown   that   the   equivalent   nodal   forces   are 
expressed  as 


yi 


z  . 

1 


N  ■  (P-  HIT  _  P   3^)d? 


"tf"n  if  +  p,  §f>« 


(B.19) 


if  the  loads  are  applied  on  an  edge  parallel  to  the 
curvilinear  coordinate  f .  The  Gaussian  numerical 
integration  is  used  to  derive  Eqs .  (B.19)  which  are  a  form 
of  Eq.  (B.7b). 

B . 6 . 2  Hygrothermal  Loadings 

The   equivalent   nodal   loadings   due   to   hygrothermal 
strains  are  given  in  matrix  form  as 


f  ■* 


y. 


[Bi]T([D2]{€°})dv 


(B.20) 


where   the   matrices   [D~]   and     {€  }     are   defined   in 
Eq.  (B.15). 
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B . 7  Element  Displacements  and  Stresses 


The  global  stiffness  matrix  as  well  as  the  equivalent 
nodal  loading  matrix  are  assembled.  Then,  the  displacements 
are  computed  by  a  frontal  solution  method.  Once  the 
displacements  are  known,  the  strains  and  the  stresses  can 
be  deduced 
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Input  the  geometric  data,  boundary  conditions, 
material  properties,  moisture  contents  and 
temperatures  at  the  nodal  points.  Generate  the 
sampling  Gauss  points. 


Calculate   all  the  element  stiffnesses  . 


Input  all  mechanical  nodal  loadings  .  Compute 
the  hygrothermal  equivalent  nodal  loadings. 


Assemble  the  equations  (B.7a).  Solve  for  the 
d  i  sp lacemen  t  s . 


Compute  the  stresses  at  all  Gauss  points. 


Fig.  B.l  Organization  of  the  F.E.M.  program. 
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o     Gauss  paint  numbers 


4    •     Element  node  numbers 


Fig.     B.2  Local     axes        f     and     n.         Gauss     point     numbers     and 

local  node  numbers  of  an  eight-point 

isoparametric    element. 
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